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We address a number of issues raised by a manuscript of Klein, Toublan, and Verbaarschot 
(hep-pli/0405180) in which the authors introduce a random matrix model for QCD with two colors, 
two flavors, and fermions in the fundamental representation. Their inclusion of temperature terms 
differs from the approach adopted in previous work on this problem (Phys. Rev. D 64, 074016 
(2001).) We demonstrate that the two approaches are related by a transformation that leaves the 
thermodynamic potential invariant and which therefore has no effect on physical observables. 

PACS numbers: 11.30. Fs, 11.30. Qc, 11.30. Rd, 12.38. Aw 

A number of model calculations of C^D at finite temperature and non-zero baryon density have revealed the 
existence of a rich phase structure A particularly interesting possibility arises when cold nuclear matter 

is brought to densities sufficiently large for quark degrees of freedom to become important. It is thought that, due 
to quark Cooper pairing, cold and dense quark matter behaves like a color superconductor Hi The system is 
characterized by a variety of color superconducting order parameters which have been studied as a function of, e.g., 
the total number of flavors, the quark masses, the quark chemical potentials, and temperature These results have 
not been confirmed satisfactorily by lattice simulations. Such calculations are difficult to perform and interpret at 
finite baryon densities. The fermion determinant appearing in the partition function is complex for non-zero baryon 
chemical potential, and standard Monte Carlo methods can no longer be applied. As a consequence, our theoretical 
understanding of the QCD phase diagram at non-zero baryon density remains incomplete. 

Random matrix theory can provide interesting elementary models that capture important aspects of the thermo¬ 
dynamic competition among different order parameters and produce robust results. Such models respect the global 
symmetries of QCD but mimic the detailed dynamics of the interactions with a form that can be treated analytically. 
When supplemented by additional physical insight regarding the form and relative importance of order parameters, 
random matrix theory can efficiently model the physics of systems with non-trivial phase diagrams. In its usual 
implementation, the fermion kinetic energy is ignored; only the interactions and their symmetries are retained. Such 
models are thus intrinsically theories of low-energy degrees of freedom which are treated at the mean field level. Near 
critical regions, random matrix results are similar to those obtained in a Landau-Ginzburg approach ill 

The effects of baryon density and temperature must be included in a manner that respects the structure of QCD. In 
the finite temperature formalism, the quark chemical potentials, pLf, and the temperature, T, enter the QCD partition 
function via the finite Matsubara frequencies tJ,f + tUn = Tf + t{2n + I)7rT (with n = 0, ±1, ±2,...). The inclusion 
of density effects through the addition of a chemical potential term, , in the Lagrangian is straightforward. 

In contrast to lattice calculations, the resulting random matrix partition function can be evaluated analytically. Its 
extremization leads to gap equations which are polynomial in the condensation fields, and which can thus be solved 
analytically or numerically. 

Different schemes have been proposed in the literature for the inclusion of temperature effects. We note first 
that, since the theory describes low-lying energy excitations, much of the critical physics is captured even if the sum 
over fermion frequencies is restricted to the two lowest terms, cUn = zLinT. When this truncation is adopted in a 
chiral random matrix theory with p./ = 0, for example, the resulting model shows a second-order phase transition 
independent of the number of flavors |^. Taking into account all terms in the frequency sum does not change this 
result fundamentally. The sole effect is the introduction of a temperature mapping that is inhomogeneous (i.e., it 
depends on T) but monotonic Hi The critical temperature is thus modified, but neither the topology of the phase 
diagram nor the order of the transition are affected. In the following, we will only consider cases where the Matsubara 
sum is truncated to the two smallest terms. 

In an early random matrix model of color superconductivity in QCD with two flavors H employed a 

partition function for which the pairing fields had opposite frequencies, w„. This choice was dictated by the physical 
observation that the pairing condensate, ('!/'^(*cr 2 )flavor(*o- 2 )coior(*o- 2 )spin t/’) (with ia 2 the antisymmetric Pauli matrix), 
is independent of time. A Fourier decomposition selects opposite frequencies for the two flavor fields. This frequency 
selection is operative independent of whether one is modeling S'(7(3)-QCD or S'C7(2)-QCD. A seemingly different 
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temperature dependence, based on the symmetries of the interactions, was recently introduced by Klein, Toublan, 
and Verbaarschot [ifl [HI ■ The temperature scheme proposed there reflects the structure of the zero-temperature 
chiral random matrix models. For SU{2) with fundamental fermions, which is characterized by real matrix elements, 
the temperature terms are chosen to be real. Similarly, both the temperature terms and the random matrix elements 
are complex for SU (3) with fermions in the fundamental representation. 

There would thus appear to be two different ways to treat temperature effects. The purpose of this paper is to 
demonstrate that the scheme based on the symmetries of the interactions is actually equivalent to that based on 
frequency selection. After summarizing the two temperature schemes in Sec. Q] we consider the transformation that 
relates them for the case of QCD with two colors in Sec. im The case of QCD with three colors is discussed in Sec. nTTi 


I. TEMPERATURE SCHEMES 


Random matrix models that are capable of studying the competition between chiral, diquark, and other condensa¬ 
tion channels are extensions of the chiral random matrix models 0 0 ll^ IT^ . In the sector of zero topological charge 
and for Nf flavors of fermions with masses m /, the partition function has the form 


Zo = n / 


Nf 




exp 






( 1 ) 


where IT is an n x n matrix which models the interaction. The elements of W are drawn on a Gaussian distribution 
of inverse variance S. Here, n is a measure of the number of low energy degrees of freedom and is to be taken to 
infinity at the end of the calculations (i.e., in the thermodynamic limit). The Dirac operator, D[IT], respects the 
chiral symmetry of QCD, {£>, 75 } = 0. In the basis of the eigenstates of 75 , D\W\ can be decomposed as 


D[W] = 




( 2 ) 


For random matrix models of QCD with SU{S) and fermions in the fundamental representation, W is complex. 
This choice corresponds to a chiral unitary ensemble characterized by a Dyson index (3 = 2. The QCD Dirac operator 
in SU{2) with fermions in the fundamental representation satisfies an additional anti-unitary symmetry, 

[C(a2)color^f,*D] =0, (3) 

with ( C (cr 2 ) color AT = 1. {C is the charge conjugation operator, 02 the antisymmetric color matrix, and K the 
complex conjugation operator.) This leads to a particular basis of states in which D is real. Accordingly, W is chosen 
real in chiral random matrix models for SU{2), and the model belongs to the chiral orthogonal ensemble with an index 
/3 = 1 H. The essential difference between the unitary and orthogonal ensembles lies in the number of independent 
random variables that are allowed per matrix element. (I.e., A complex number has two degrees of freedom; a real 
number has only one.) This difference is important in determining the statistical properties of the Dirac operator 
and, in turn, influences the phase diagram. 

The inclusion of quark chemical potentials, /r/, and the teimerature, T, follows from the frequency representation 
of the partition function in the finite temperature formalism [^ ll fil | . The path integral is evaluated for fermion helds 
that obey antiperiodic boundary conditions in Euclidean time. Thus, the Fourier decomposition 

V'('r) e"“"'"V'(wn) (4) 


contains only the odd Matsubara frequencies 

uin = i{2n + l)nT. 

The dependence on fj,f and T appears as a component of the Lagrangian of the form 

Li= •!/'/74(9r + M/)V’/, 

/ 

whose Fourier representation is given as 



( 5 ) 

( 6 ) 


( 7 ) 
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There is a certain freedom associated with the introduction of such a term in random matrix models. Equation 0 
suggests that we should add a deterministic term with the block structure of 74 to the Dirac operator of Eq. ©■ 
With Euclidean matrices, this gives 


D[W]^ 


0 

iW^ +LO-{^lf,T) 


lW + Lo+{^Jif,T)\ 

0 ) 


( 8 ) 


where are n x n matrices. Here, n is a compound index which combines random matrix and Matsubara frequency 
indices. Because the ensemble of matrices W is invariant under unitary (/3 = 2) or orthogonal (/3 = 1) basis 
transformations, many equivalent forms of lo± can be adopted, and their detailed form cannot be specified uniquely. 
(We will also consider another transformation below.) We know, however, that the eigenvalues of ui± must be those of 
the operator in Eq. 0 . In the following, we restrict ourselves to the two lowest Matsubara frequencies, ujn = iinT. 
The appropriate forms of uj± must then have two series of n/2-degenerate eigenvalues ± zttT. 

We now turn to the various schemes proposed for implementing temperature dependence. In an early model of 
diquark condensation |3) 0 in QCD with two flavors, we considered random matrix interactions that included spin 
and color quantum numbers explicitly. This choice leads to an expanded block structure 


3 3 

kE = ^ ^ (cr^ (g) (Ta) Af,a, 

^—0a—1 


(9) 


where cr^ = (1, ia) represents spin degrees of freedom with the Pauli matrices a, while da are the SU{2) color matrices. 
The interactions are to be integrated over the real matrices whose elements are drawn on a Gaussian distribution. 
The W appearing in Eq. © is complex independent of the number of colors. Because the dependence on spin and 
color is deterministic, the model nevertheless has the appropriate number of degrees of freedom per matrix element. 
For example, for Nc = 2 the Dirac operator D[IT] in Eq. |(21) with the block structure of Eq. satisfies the antiunitary 
symmetry of Eq. ©■ More generally, we have shown that for = 3 and Nc = 2 the eigenvalues of the Dirac operator 
ZI[IE] exhibit the spectral correlations expected for the chiral unitary and orthogonal ensembles, respectively [T3 |. 
The deterministic part of the Dirac operator is given as 


UJ+{fif,T)=UJ-{nf,T)= (Ispin (g Icolor) (m/ + , (10) 

where the choice of s/ = ±1 follows a selection rule. First, s/ is taken positive for flavor / = 1 and negative for flavor 
/ = 2. The contribution from the opposite choice, si = —1 and S 2 = 1, is then added to obtain the thermodynamic 
potential. (In Refs. 0,0, we considered the limit fj,i = fj. 2 ] in that case, the two contributions are identical and the 
sum over pairs amounts to a factor of two.) This procedure is suggested by the Fourier expansion of the condensates 
in a microscopic theory. The chiral condensate is given by the equal time correlator Cu ~ {tjj{T)'4’{T)), which is also 
time independent. A frequency decomposition of the fields, Eq. 0 , leads to an expression which contains a sum over 
all frequency pairs. Since the correlator Ca is time independent, the only non-vanishing terms are those for which 
the Fourier components ip{uJn) and ip(iiJm) have equal frequencies (or ra = n): Ca T{J2n4’{^n)''P{^n))■ Similarly, 
the diquark condensate is given by the equal-time correlator Ca ~ (V'^(''')(*0’2)coior(*0-2)spin(*0-2)flavor V'C''")) s-iid is 
constant in time r. A Fourier decomposition now selects opposite frequencies for the two fields, which necessarily 
correspond to different flavors because of (*cr 2 )flavor- This mechanism was implemented in the random matrix model 
through the selection of appropriate pairs (si, 52 )- 

In recent work mini, Klein, Toublan, and Verbaarschot treated the chiral random matrix ensembles, studying 
diquark condensation for the case Nc = 2 only. The zero temperature and zero fij version of their model follows 
Eq. 0 with W complex for the unitary ensemble and real for the orthogonal ensemble. In their approach, the 
dependence on fj,f and T follows the symmetries of the chiral ensembles. For the unitary ensembles, lo± are complex. 
In the appropriate sub-basis of frequency states, they assume the diagonal form 

= (“) 

where each block matrix has dimension n/2 x n/2. For the chiral orthogonal ensembles, the antiunitary symmetry 
of Eq. 0 is imposed, which leads to a particular basis for which the ui± matrices are real. Ensuring that their 
eigenvalues are given as ± inT, the proposed dependence is 


UJ+ 


( A*/ -^T 

\ Ilf 


UJ- 


( 12 ) 
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where the dimensions of the block matrices are again n/2 x n/2. 

At first sight, the two forms of uj±, Eqs. HU and eg, might seem non-equivalent and also different from that 
of Eq. HI d|l . There thus appears to be a disagreement between the temperature dependence that is obtained from a 
generic structure suggested by symmetry and one which is based on physical knowledge of the frequency couplings 
that are allowed. In the next section, we consider the case Nc = 2 and show that the dependence expressed by the 
selection rule in Eq. m leads to the same effective thermodynamic potential as that obtained with the dependence 
of Eq. (I12II . We comment on the correct form to be used for Nc = 3 in Sec. IIIII 


II. EQUIVALENCES BETWEEN THE TEMPERATURE SCHEMES 


It is useful to compare the temperature schemes of Eqs. (nni, mi, and (EJ in a common framework. To this end, 
we recast the model developed in [^||| in a form close to that of Refs. Eai [lli| . We consider a model for Nc = 2, two 
flavors, and zero quark masses. We work with an interaction W having the expanded block-structure of Eq. and 
derive the thermodynamic potential for the partition function of Eq. O- After integrating over the random matrices 
Afj_a and performing a Fierz transformation, we obtain the partition function 


^ = n /exp (^) exp j , (13) 

where Y represents the four-point interaction. If only the chiral and diquark channels are retained, 

^ (V’li • V'li) + ■ ^’^ 2 ) (V'l2 ■ i’L2) + 

(V’L • Pa ■ (V'l 2 • Pa ■ V'li) + (V-L ' (V'L ' Pa ■ V-m)) ■ (14) 

Here, we have adopted the notation (j) ■ x = X)r=i 4’iXii where i is the random matrix index, and have introduced the 
projector Pa = (* cr 2 )^n(* 0 - 2 )color- 

Following and pTl| . bosonic auxiliary variables are introduced via a Hubbard-Stratonovich transformation that 
accounts for the different combinations of spin, color, and flavor quantum numbers. This transformation is given by 

(15) 


where A is a complex, antisymmetric 16 x 16 matrix, and y{A) is the bilinear 



(16) 


Here, 1„ represents the identity matrix in the combined space of random matrix/frequency indices and B is given as 


B = 



In a basis of states = {'ipi , 1 / 12 ) V'l j ^' 2)1 Ih® matrix A assumes the form 

/ 0 ArPa at 0 \ 

_ I -ArPa 0 0 a* \ 

-at 0 0 -AIPa ’ 

V 0 -a* AIPa 0 / 


(17) 


(18) 


where each block is 4 x 4 (two spins x two colors), and the auxiliary variables A^, Ar, cti, and CT 2 are complex scalars. 
The temperature and chemical potential terms of Eq. m can be arranged as 




0 UJ+ 

ui- 0 





( 19 ) 
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Combining Eqs. Ca), 113), and 113 and integrating over the fermion fields, we finally obtain 



For later discussions, it is important to note that the diagonal blocks of the argument of the logarithm are given as 
an exterior product with 1„. Different Matsubara frequencies will thus correspond to the same auxiliary matrix, 
or B. 

Equation is the central expression which we will use to compare the different temperature schemes. We will 
concentrate on the same condensation channels as those considered earlier 1^. We take equal quark chemical 
potentials, fJ-i = fJ ,2 = The chiral condensate is considered to be flavor independent so that ai = a 2 = cr. Moreover, 
we have {'^'^rPaiPir) = and the auxiliary fields satisfy the equality Ar = Ar = A. 


A. Selection rule scheme 


The temperature dependence proposed in Refs. Hi can now be represented as (/ii = /r 2 = /r) 

/ /r + iirT 0 \ 

\ 0 u — iirT /,, , , , 

> ■ rr. r. \ Matsubara (21) 

I fJ- — *7ri 0 \ 

V 0 ^ + Matsubara' 

Inserting this form in Eq. leaves us with a straightforward evaluation of a trace. We obtain 

Z = j VA (22) 


f UJ+{^il,T) = w_(^i,T) 

I w+(^2,T) = W_(^2,T) 


with an effective thermodynamic potential 

D(cr,A) = (cr^ + I Ap) - log ((ct -/r)^ + I Ap+ 7r^T^) - log ((cr + ^)^ + I Ap TTr^r^) . (23) 

(The inclusion of a quark mass mf = m would amount to replace cr by cr + m in the argument of the logarithm.) In 
the thermodynamic limit n —>■ oo, one can evaluate the partition function integral exactly by a saddle point method 
and obtain 


^ g 4nQ((T3j^udle lAsaddle) 

n—^oo 


(24) 


where the fields crgaddie and Agaddie are solutions of the gap equations 

0, (25) 

0 . (26) 

Although the potential in Eq. (EHll has been obtained using a formalism different from that used in Ref. |3 , we recover 
the earlier result. We now turn to other temperature schemes. 


dn 

on 


(ct,A) = 
(ct,A) = 


B. A selection rule scheme with explicit flavor symmetry 


At first sight, the temperature dependence in Eq. might seem to violate flavor symmetry. However, this is not 
the case. The sign reversal in T for flavor 2 can also be expressed as a permutation of the two Matsubara frequencies, 
since 




iirT 0 
0 —inT 


iirT 0 
0 —iirT 


(27) 
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Therefore, we can go to an equivalent representation by permuting the frequency states for flavor 2. The selection of 
opposite frequencies applies to diquark and antidiquark bilinears for which the projector P/^ is sandwiched between 
the flavor fields as in Eq. da- Hence, the effect of permuting positive and negative frequencies for flavor 2 transforms 
Pa into 

Pa ^ Pa — Pa f ^ Q j — (^cr2)spin ^ (^^2)color ^ f Q j 

\ / Matsubara \ / Matsubara 

By contrast, the chiral bilinears are unaffected by the frequency permutation, e.g., = '0+^^'+ + The 

sole effect of permuting the frequencies for flavor 2 is thus to make the frequency selection mechanism appear in 
the structure of the order parameter through the projector Pa- The frequency component of this operator is flavor 
diagonal and Pa therefore respects the required flavor symmetries. 


C. Antisymmetric temperature dependence 


The apparent disagreement between a temperature dependence based on symmetry and one based on selecting 
frequency couplings would immediately be resolved by a transformation relating Eqs. 03 and m- In the antisym¬ 
metric representation of Eq. 03 , the basis of states is meant to be that for which all matrix elements of W are 
real. This is certainly not the basis set to be used in our model, Eq. 03 , since the spin and color dependence of 
Eq. @ makes W complex. Nevertheless, the Dirac operator P[IE] of Eq. (|3) satisfies the anti-unitan^ symmetry 
of Eq. (|3J), and we could in principle choose the appropriate spin and color basis to obtain a real W Q- This is a 
consequence of the fact that the matrix elements contain the proper number of random degrees of freedom. In the 
following we are interested in the transformation properties of the deterministic temperature terms u;± and thus do 
not consider transformations in spin/color space. We focus instead on the effects of a partial transformation acting 
only on frequency indices. 

We wish to establish the relationship between the temperature dependence based on the selection rule of Eq. 03 
and an antisymmetric dependence. We start from the representation of Sec. cm for which 

w± = M + ( Q _^ , (29) 


and the diquark order parameter contains the projector Pa of Eq. (I28II . We rewrite the diagonal temperature form as 


f inT 0 \ 

0 -inT) 


= Urr 


0 nT 
-nT 0 


uL 


(30) 


where 


^"^“72 ( 

vl * ) 

(31) 

is unitary. We then perform the transformation 



= 


(32) 

V'fl.L = 


(33) 

with 



Tj f Ujn 0 ln/2 0 

^ 0 Ispin ^ Icolor ^ Iflavor- 

(34) 


This transformation leaves the partition function of Eq. 03 unaffected. Since U only acts in frequency space, it does 
not affect the quadratic term, TrAA0 The logarithmic term is also unaltered since TrlogM = logdetM, and the 
determinant is invariant under the transformation 


M ^ M 


0 

0 W 


M 


U 0 \ 

0 U* )■ 


(35) 


Hence, the transformation, U, brings the temperature representation of Sec. IIIBI to an equivalent representation with 
an antisymmetric temperature dependence. 
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We now examine the transformation of the bilinear form of Eq. oa bearing in mind that we started with the 
representation of Sec. ITTBl for which the diquark terms contain the projector of Eq- (EEJ- Mesonic bilinears 
of the form V’lV'j (where i and j represent some combination of spin/color/flavor indices) are invariant under the 


P-transformation as 



(36) 

By contrast, diquark bilinears of the form P^'i’j s-re transformed to 


iff pA'lp 3 = ip^PApy 

(37) 

with 


p'a = UIPaU^ = Pa, 

(38) 


where again, Pa = (*cr 2 )spin(*cr 2 )coior is block diagonal in frequency space. 

To summarize, the temperature dependence in the new basis is antisymmetric, while the diquark order parameter 
no longer contains an explicit frequency selection operator. This is precisely the form proposed in Ref. EH- We thus 
have shown that it is equivalent to the selection rule of Eq. and, by extension, to the one originally proposed in 
Ref. HH. To complete the proof, it is readily verified that the various temperature schemes of Eqs. (fTn)l. jEJ, and 
(P|l . used with the corresponding forms for the diquark projector (Pa, Pa, or P^ = Pa), all yield the thermodynamic 
potential of Eq. 


D. Pion condensates 

Although the temperature dependence in Eq. m produces the correct result for modelling diquark and chiral 
condensation, it is inconvenient if one also wishes to consider the possibility of forming a pion condensate, ('0^75 (jcr2)'*/')- 
Pion bilinears must be taken between fields of equal frequencies and opposite flavors, whereas elements associated 
with diquark bilinears require opposite frequencies and opposite flavors. Considering both channels thus leads to 
some tedious bookkeeping. However, the equivalent form of Eq. (EHli can be implemented much more easily and 
extended to the inclusion of pion bilinears. The resulting form takes a diagonal temperature dependence, considers 
diquark bilinears with an off-diagonal projector Pa, while it leaves both the chiral and pion bilinears diagonal in 
frequency space. We now show that such prescription gives the same results as those obtained with an antisymmetric 
temperature dependence. 

A Fierz projection of the random matrix interactions on the pion condensation channel yields a four fermion 
component 

Ett = ('*/'^75(*O'2)flavor'0) . (39) 

Using again a Hubbard-Stratonovich transformation to introduce an associated auxiliary variable, p, one is led to 
consider the auxiliary matrix 


/ 0 A/jPa a* P \ 

-ArPa 0 -p 0-2 I 

-at p 0 -AIPa ■ 

V -p -at A*Pa 0 / 


(40) 


As previously, the frequency selection rule is introduced by replacing Pa by Pa in the diquark and antidiquark 
bilinears of Eq. OS or equivalently in the argument of the logarithm of Eq. (HOI). Such a temperature scheme is 
completely consistent with the spirit with which we derived the temperature dependence in Refs. |^; we have 
verified that it produces the same thermodynamical potential as that obtained with an antisymmetric temperature 
dependence. The result is identical to that of Ref. pTl| . 


III. DISCUSSION AND CONCLUSIONS 

For Nc = 2, we have resolved the apparent disagreement between the dependence based on an antisymmetric 
temperature term and that based on a frequency selection rule. The two representations are related by a unitary 




transformation which leaves the thermodynamic potential unchanged and which therefore does not affect physical 
observables. The change of basis is, however, reflected in the structure of the diquark order parameter or, equivalently, 
in the expression for the frequency selection rule. 

Technically, temperature enters through the imposition of an antiperiodic boundary condition on the fermion fields. 
There is a certain degree of freedom in representing temperature terms in a random matrix model. However, as soon 
as one imposes a given block structure (e.g. diagonal or antisymmetric), one is actually selecting a specific basis for 
the frequency states. As a result, attention must be paid to the physical meaning of the basis states when solving 
for the partition function and, accordingly, in the structure of the order parameter. If a change of basis appears 
useful for mathematical convenience, one must naturally keep track of the effects of such transformation on the states 
and transform the order parameter accordingly. From a general point of view, any transformation with the block 
diagonal form of Eq. m with Um unitary leaves the thermodynamic potential invariant. There are thus many 
equivalent Matsubara bases which can be chosen and correspondingly many correct associated forms of uj± and of 
the diquark order parameter. Among these forms, the antisymmetric one appears to be particularly convenient since 
the corresponding diquark projector is block diagonal in frequency space. The price to be paid for working with a 
different form is merely the need to rely on the knowledge of the physical meaning of the basis states and on which 
frequency couplings are allowed. This will determine the form of the corresponding diquark projector. 

The situation is not very different for QCD with = 3. If one chooses a diagonal temperature term, one still 
needs to select opposite frequencies for the fields involved in the diquark condensate. In this case, however, the 
frequency selection mechanism cannot be revealed by symmetry arguments alone since QCD with N^. = 3 does not 
satisfy an antiunitary symmetry. As above, there are many correct forms of io± which can be used together with the 
corresponding prescriptions for the frequency couplings. In Ref. p7i|. the authors used the diagonal form of Eq. m 
but did not consider diquark condensates. Thus, there was no need for specifying a diquark projector, and the choice 
of a diagonal form of lo± was as correct as any other. On the other hand, diquark condensation was investigated in 
Ref. using a diagonal form. There, it was necessary to supplement the temperature dependence with an explicit 
frequency selection rule. 

In conclusion, the proper temperature dependence to be used in a given random matrix model depends primarily 
on the form of the condensates that may develop. Symmetry arguments are useful in identifying the mechanism for 
selecting appropriate frequencies, but they are not always sufficient. 
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